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Abstract 

. 12 Ruthenium-based perovskite systems are attractive because their structural, electronic and mag- 

13 netic properties can be systematically engineered. SrRu03/SrTi03 superlattice, with its period 

o 

14 consisting of one unit cell each, is very sensitive to strain change. Our first-principles simulations 
. 15 reveal that in the high tensile strain region, it transits from a ferromagnetic (FM) metal to an 

>: 

, i6 antiferromagnetic (AFM) insulator with clear tilted octahedra, while in the low strain region, it is 

ON 1 17 a ferromagnetic metal without octahedra tilting. Detailed analyses of three spin-down Ru-t2g or- 

in; 

qq ■ is bitals just below the Fermi level reveal that the splitting of these orbitals underlies these dramatic 

o; 

£N| . 19 phase transitions, with the rotational force constant of RuOg octahedron high up to 16 meV/Deg 2 , 

t-H ■ 

20 4 times larger than that of TiOg. Differently from nearly all the previous studies, these transitions 

• i— i ■ 

^ ■ 21 can be probed optically through the diagonal and off-diagonal dielectric tensor elements. For one 

5r ■ 22 percent change in strain, our experimental spin moment change is —0.14 ± 0.06 ub, quantitatively 

23 consistent with our theoretical value of —0.1 ub- 

24 PACS numbers: 74.70.Pq 68.65. Cd 75.70.Cn 78.67.Pt 
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Strontium ruthenate [l| belongs to a broad scope of perovskites and has attracted ex- 
tensive attention, due to its exotic properties. When slab thickness decreases, its itinerant 
ferromagnetic (FM) phase disappears Qj. This has motivated many experimental 3^ 6 ] and 
theoretical investigations [?-9] to search the origin of the loss of the FM metallic phase. 
Another perovskite, SrTiC>3 (STO), is a good insulator, and also has attractive physical 
properties such as superconductivity and two-dimensional electron gas on its surface and 



interfaces 



id, 



It would be a fascinating idea to investigate a superlattice which consists 



of one layer of "conducting" SrRuOs and one "insulating" layer of SrTiOs, or SRO/STO su- 
periattice. E„tall y , thick—lattice, have already heen fabricated QQ. The 



10 Curie temperature decreases 



14 



Il5j | with the decrease of the period of the superlattice. The 



magnetic moment of the Ru atom is suppressed, and no FM ordering was identified in a 
1/1 superlattice. A question is thus raised whether there is indeed any magnetic ordering 
established in the SRO/STO 1/1 superlattice. 

With the advent of the state-of-the-art molecular-beam epitaxy (MBE), pulse laser de- 
position (PLD) and other growth techniques, it is now possible to fabricate interfaces with 
atomic sharpness. This is particularly true for systems with a very small lattice mismatch, 
such as SRO/STO. More importantly strain can be controlled through different substrates. 
For instance, Ti02 and MgO 16[ can induce an in-plane strain of -4% and 7%, respectively. 



Piezoelectric substrates like PMN-PT even allow one to control the strain in real time 
This greatly facilitates materials engineering. 

In this Letter, we report that such a superlattice undergoes an intriguing magnetic 
phase transition under epitaxial strain. Our experiment first shows that in a well prepared 
SRO/STO superlattice sample, ferromagnetic ordering can survive down to ultrashort period 
(1/1). We can tune its magnetic properties by applying different strains. Our first-principles 
calculations further reveal that the in-plane strain can drive the system from a ferromag- 
netic to an antiferromagnetic phase at a critical strain £ c = 5%. Within the ferromagnetic 
phase, three structurally different phases are identified: below 0.25% (a phase), the RuOe 
and TiOe octahedra rotate in the opposite direction but without tilting; between 0.25% and 
2.5% ((3 phase) tilting starts and the rotation angles of both RuO@ and TiOg are reduced; 
and above 2.5% (7 phase), both RuOg and TiOg rotate in the same direction. To understand 
these dramatic changes, we carefully examine the borderline between phases and find that 
the frontier spin-down Ru-t2 g orbital is mainly responsible for the phase transition, where 
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1 its occupation changes with the strain. These phase transitions directly lead to a qualitative 

2 difference in dielectric tensor, a signature that can be probed experimentally. 

3 We grew [SRO/STO]3o superlattice samples using laser MBE on three different sub- 

4 strates: STO, Nb:STO and LaSrAlTaO (LSAT). The growth of 1/1 superlattice films was 

5 monitored by reflection high energy electron diffraction (RHEED), which showed a layer- by- 
e layer growth mode, (see the Supplementary Materials for more details.) Lattice mismatches 
? in terms of STO lattice constant for these substrates are 0%, 0.05% and -0.92%, respec- 
s tively. The structures were characterized by x-ray diffraction (XRD) using the synchrotron 



9 radiation [18| beamline BL14B1 of Shanghai Synchrotron Radiation Facility (SSRF), shown 

10 in Fig. 1(a). The data shows unambiguously that the superlattices are smooth and free 

11 of any second phase. From the (002) peaks (see Fig. (Ha)), we estimate the out-of-plane 

12 lattice constants for these films to be 3.982 A, 3.978 A and 4.009 A, respectively, meaning 

13 that the samples are strained according to the substrate lattice. The Laue oscillations of 

14 the peak due to the total crystalline film thickness indicate a good [001] orientation of the 

15 film. Surface atomic force microscopy (the inset of Fig. (Ha)) reveals that all these samples 

16 have a smooth termination with the roughness below 1 uc. 

17 The magnetic properties were measured at 5K by vibrating sample magnetometer (PPMS 
is VSM Option Release 1.2.4 Build 1). Although the magnetization is suppressed compared 

19 to the bulk SRO, clear hysteresis loops (M-H) are observed in all our samples. Superlattice 

20 recovers some of the lost FM ordering in SRO thin films. To our knowledge, this has not been 

21 reported before at 1 uc-thick SRO thin films and superlattices. Though the coercive field 

22 weakly depends on substrates, both the remanence and saturation field strongly depend on 

23 the substrates. The spontaneous magnetization changes from 0.12 ± 0.03 hb (on Nb:STO) 

24 to 0.30 ± 0.03 hb (on LSAT) per Ru. The Ru magnetic moment decreases as the strain 

25 increases (see the inset in Fig. 1(b)). For one percent change in strain, the magnetic moment 

26 change is -0.14 ± 0.06 Since film thickness and growing conditions are the same for all 

27 our samples, the lattice mismatch, or the strain, is directly responsible for the magnetic 

28 properties change. 

29 To understand the strain effects in this superlattice, we resort to first-principles calcu- 

30 lations. For bulk SRO, extensive calculations have been performed to investigate changes 

31 in structural, electronic, and magnetic properties {§], Il9l-22|. but very few on a superlattice 
23j-|25|. We carry out first-principles calculations on a STO/SRO superlattice (see the inset 
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1 in Fig. 1(a)) within t 

2 repulsion (LSDA+U) 



le local spin density approximation plus Hubbard on-site Coulomb 



18 



261 ] . Within this scheme, the magnitude of U e g is treated as an 



3 empirical parameter, which will be discussed in the following text and the Supplementary 

4 Materials. 

5 The in-plane tensile strain £, defined as £ = (a— ctsTo)/ a STO; is applied to the superlattice. 
e Here a is the lattice constant in use while «sto is that of the parent compound. The in-plane 
7 strain is changed from -4.5% to 6%. At each strain, the ion positions and the out-of-plane 
s lattice constant are fully relaxed. There are three distinctive angles to characterize the 

9 structure change. One is the tilting angle <f> around an axis parallel to (001) plane (see the 

10 inset in Fig. QJc)). The other two angles #r u and 6^ denote rotations around [001] axes of 

11 the RuOg and TiO e octahedra, respectively. 

12 Figure [T](c) shows the energy difference between FM and AFM phases as a function of 

13 strain. In a wide range of strain, the FM alignment between the neighboring Ru atoms 

14 is energetically favored. The magnetic moment is 2 \xb per SRO formula unit, which is 

15 consistent with the calculated ferromagnetic ground state of bulk SRO 22j . However, when 

16 the tensile strain exceeds ~5%, the structure with c-type AFM phase becomes more stable 

17 than that of FM. Therefore, a magnetic phase transition occurs at this critical strain. One 
is notices that the total energy difference between the AFM and FM structures is linear with 

19 respect to strain but with three different slopes (see the caption of Fig. 1). This suggests 

20 structural phase transitions taking place as a consequence of strain variation. 

21 The observation is indeed verified. Structurally, the FM region can be subdivided into 

22 three different phases: £ < 0.25% (a phase), 0.25% < £ < 2.5% {(3 phase) and 2.5% < 

23 £ ^ 5% (7 phase) (see Figs. 1(c) and 1(d)). Within a phase the neighboring octahedra 

24 rotate counterclockwise with respect to each other, where is always positive but 6ti is 

25 negative (see Fig. 1(d)). There is no tilting, i.e. = 0°. The quenching of tilting renders 

26 the superlattice with a high symmetry of P4/mbm, making it easy be detectedjsee below). 

27 Similar observations have been reported experimentally in single crystals (si, 0-29]. Once 

28 the strain exceeds 0.25%, the symmetry is reduced to P2\/c. The phase features a tilting 

29 {(f) 7^ 0°) and two octahedra TiOg and RuOg rotating in the opposite direction (see 6^ < 

30 and #R U > 0). In the 7 phase the TiOq octahedron rotates in the same direction as RuO@. 

31 To shed light on this structural phase transition, we zoom in a small region around 

32 £ = 0.25% (see the small dashed box in Fig. [0(d)). We manually change four structural 
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parameters (0, 6*r u , #ti and £) around their respective equilibrium values while keeping the 
rest unchanged. To make a quantitative comparison, we choose the structure at £ = 0.45% 
with a tilting angle = 3.88° as the reference structure since it is near the critical point. 
The energy difference curves are plotted in Fig. 12(a) . The energy minimum is indeed at 
its global minimum, since each curve has a minimum at 0°, with a small deviation in A# Ti 
due to our energy threshold of ±1 meV in our optimization procedure. Contributions from 
each angle are very different. The rotational angle A6*r u has the strongest effect on the 
energy change among all the angles. If we fit the energy change to a harmonic potential 
|A"a0 Ru (A#r. u ) 2 , we find K^g Ru = 16 meV/Deg 2 , or 53 eV/rad 2 [30j. This is much larger than 
that of Ka</> = 1 meV/Deg 2 and K&g Ti = 4 meV/Deg 2 . In addition, the energy change of 
A6*Ti is highly anharmonic. We expect experimentally Raman spectra can easily distinguish 
them. 

It is conceivable that the above strongest contribution from A#r u must be associated with 
the electronic structure of Ru ions. To see this, we integrate the Ru-t2 g projected density of 
states from -10 eV to the Fermi level. Figure 2(b) shows that the occupation of these orbitals 
changes substantially with strain. If we compare its change with the phase change in Figs. 
1(c) and 1(d), we find a very nice match between them. Starting from the strain away from 
the critical one, the occupations in d zx and d yz are similar. Near the phase boundary of the 
structural phase transition their occupations differ from each other. The d yz orbital gains 
electrons while the d zx loses electrons. After tilting sets in, the occupations of these two 
orbitals again become close to each other. No other element has these characteristic changes. 
Moreover, since the Fermi level is mostly contributed by the Ru-0 hybridized states, this 
explains why the energy change with respect to the rotation of RuOg is most pronounced. 
At the FM/AFM phase transition point, electrons are transferred from the d yz / zx to the d xy 
orbital. One spin down electron resides almost entirely in the d xy orbital. More details are 
provided in the Supplementary Materials. 

We also investigate the effect of the Hubbard U on the above phase transitions. Struc- 
turally, U e fi has a minor effect 8j. Changes in the rotation angles at strain £ = 1% and 
£ = 4.5% are too small to show. The largest change of about 2 degrees in the RuOg rotation 
(Fig.2(c)) is found at a highly strained case (£ = 6%). But none of these affects the above 
phase separation. Therefore, the PA/mbm P2i/c transition is robust [3]. Figure 2(d) 
shows that the band gap for the AFM phase is established for U e s > 1 eV, and the total 
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i energy favors the AFM phase for U e g > 1.5 eV. This is expected since it is well known that 



3l|. 



2 the strong on-site correlation favors an AFM phase 

3 Since U c g effectively splits and shifts band states, the spin moment is fixed. To compare 

4 with our experimental spin moment change, we set U e g to zero [32| and compute the spin 

5 moment change. The results are shown in the inset of Fig. WO 3 )- For every percent change 
e in strain, our theoretical magnetic moment change of -0.1 fiB agrees with our experimental 
7 value of -0.14 ± 0.06 \ib quantitatively. 

s Finally, we demonstrate that both predicted structural and magnetic transitions are de- 

9 tectable optically. The structural phase transition at £ = 0.25% breaks the mirror symmetry 

10 (C2v), while the magnetic ordering transition changes the band structure. The former leads 

11 to a dramatic difference in the off-diagonal element of the dielectric tensor, and the latter 

12 leads to another big difference in the diagonal elements. In other words, we probe structural 

13 and magnetic phase transition using two different tensor elements. Since VASP does not 

14 include the intraband contribution, we decide to use Wien2K to compute the tensor since 

15 both interband and intraband transitions are taken into account 18]. 

16 The off-diagonal element of the dielectric tensor is used to probe the first order structural 

17 phase transition. In a phase, C% v mirror symmetry exists and all the off-diagonal elements are 
is zero. When tilting sets in and the mirror symmetry is broken in /3 phase, nonzero e yz emerges, 

19 a manifestation of the beginning of octahedra tilting (see Fig. |3]^a) and more discussions in 

20 Supplementary Materials). At the critical point the metal-to-insulator transition occurs. To 

21 observe this transition, the diagonal element is used. In the metallic phase, the low energy 

22 excitation consists of a plasma contribution, and the diagonal elements of dielectric tensor 

23 diverge like 1/u. In the insulating phase, on the other hand, only interband transitions are 

24 left. The low energy divergence no longer exists. Such change can be seen in Fig. [3]^b). 

25 Thus the two major phase transitions can be detected through a simple optical setup. 

26 In conclusion, we have shown that SRO/STO superlattice preserves its ferromagnetic 

27 ground state at ultra-short limit. Our experiment has demonstrated that their magnetic 
2B properties are tunable through different strains induced by different substrates. This is 

29 confirmed in our theory. Our theory further reveals a strain-dependent phase evolution for 

30 SRO/STO superlattice, where increase in strain can drive the superlattice from a ferromag- 

31 netic metallic phase to an antiferromagnetic insulating phase. There are three phases within 

32 FM. In the a phase, the RuOg and TiOg octahedra do not tilt, but in the /3 and 7 phases, 
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they do. We have shown that the Ru-t2g orbital underlies these multiple-facet changes, which 
can be detected experimentally. By examining the effects of Hubbard U, we find that our 
theory with LSDA+U does qualitatively support our experimental findings that the strain 
induces changes in magnetic properties. Therefore, our findings are significant as they reveal 
fascinating opportunities in the Ru-based strongly correlated electronic systems, which are 



crucial for future applications in ferroics and nano devices [33|-|35|. 
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FIG. 1: (Color online) (a) X-ray diffraction pattern for the sample grown on STO substrate. Left 
inset: Superlattice structure, where Sr, Ru, Ti and O atoms are shown in green, brown, blue 
and red, respectively. Right inset: Surface AFM image for this sample, (b) Hysteresis loops for 
[SRO/STO]3o superlattices grown on STO, Nb:STO, and LSAT substrates, respectively. The inset 
shows the magnetic moment of the Ru atom as a function of strain. (Exp: unfilled symbols and 
theory: filled symbols.) (c) Theoretical in-plane strain dependence of the total energy difference 
between AFM and FM phases. Three phases a, {3 and 7 in FM phase are highlighted by shaded 
in three different colors. The slopes for the three fitting lines are -0.03, -0.01, -0.04 eV/percent 
in strain, respectively. Inset shows the tilting angle (f> and two rotational angles #r u and #Ti- (d) 
Optimized 4>, #r u and #Ti as a function of £. The dashed box is the region that is further examined 
in Fig. 2(a). 
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FIG. 2: (Color online) (a) Relative total energy as a function of the rotational and tilting angles. 
Here A0 = (f> - 4>° , A0 Ru = 9 Ru - 6> Ru , A6> T i = #Ti - #Ti> where 0°, 6> Ru and 0^ are their respective 
equilibrium values at the critical point 0° = O°,0 Ru = 10.9°, 0!^ = -4.7°). E is the energy of 
the equilibrium structure. A0Ti has its lowest point at +0.6°, but this is within our relaxation 
accuracy of 2 meV. (b) Spin-down occupation of the three Ru t2 g orbitals. See spin-up occupation 
in the Supplementary Materials. The occupation is calculated by integrating the projected DOS 
from -10 eV to the Fermi level within the Wigner-Seitz radius of 1.402 A. (c) Tilting and rotation 
angles at £ = 6% as a function of U e g. (d) Band gap for the AFM phase and total energy difference 
as a function of U c ff. 
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FIG. 3: Different optical tensor elements are used to monitor structural phase transition (from an 
untilted to a tilted structure), and the metal-to-insulator transition (from a metallic FM phase to 
an insulating AFM phase), (a) Off-diagonal element, Im(e yz ), as a function of photon energy Hlo. 
The untilted structure has a null signal, while the tilted one has a signal, (b) Diagonal elements, 
Im(e xx ), as a function of Huj. The focus is on the lower energy side. In the metallic phase the 
Im(e xx ) diverges, while no divergence exists in the insulating phase. 



12 (August 30, 2012) 



